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ABSTRACT 

We report results from the Wyoming Survey for Ha (WySH), a comprehen- 
sive four-square degree survey to probe the evolution of star-forming galaxies 
over the latter half of the age of the Universe. We have supplemented the 
Ha data from WySH with infrared data from the Spitzer Wide-area Infrared 
Extragalactic (SWIRE) Survey and ultraviolet data from the Galaxy Evolu- 
tion Explorer (GALEX) Deep Imaging Survey. This dataset provides a multi- 
wavelength look at the evolution of the attenuation by dust, and here we com- 
pare a traditional measure of dust attenuation {L{T I R) / L{FUV)) to a diagnos- 
tic based on a recently-developed robust star formation rate (SFR) indicator, 
[Haobs + 24/xm]/Haobs- With such data over multiple epochs, the evolution in 
the attenuation by dust with redshift can be assessed. We present results from 
the ELAIS-Nl and Lockman Hole regions at ^ ~ 0.16, 0.24, 0.32 and 0.40. 
While the ensemble averages of both diagnostics are relatively constant from 
epoch to epoch, each epoch individually exhibits a larger attenuation by dust for 
higher star formation rates. Hence, an epoch to epoch comparison at a fixed star 
formation rate suggests a mild decrease in dust attenuation with redshift. 
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INTRODUCTION 



A fundamental property of a galaxy is the number of stars it gives birth to each year, 
commonly known as its star formation rate (SFR). Likewise, the average number of stars 
forming per year per cubic megaparsec, the cosmic star formation rate density (SFRD), is a 
key parameter in the study of galaxy evolution. To quantify these fundamental parameters, 
a number of different star formation rate indicators are used, and they stem from a wide 
range of wavelengths. The star formation rate can be estimated using the X-ray, ultraviolet 
(UV), infrared (IR) or radio continuua, a var i ety of Hydrogen recombination lines, or heavier- 
atomic forbidd en lines (e.g.. lKennicuttlll998l : ICardiel et al.ll2003Hlglesias-Paramo et al.ll2004j : 
Ly et al.ll2007l :and references therein). Three of the most widely-used star-forming measures 
are the ultraviolet and infrared continua and the Ha recombination line, which respectively 
probe i) the emission from hot, relatively young (OB) stars, ii) starlight that has been re- 
processed by interstellar dust and iii) the H II regions. Because these three star formation 
rate indicators arise from three different processes and e mit over three different wavelength 



regimes, there are know n discrepancies between them (IHopkins et al.l l2001at iTresse et al. 



2OO2I : iTeplitz et al.ll2003l ). The observed star formation rates will clearly disagree due to the 
variable effect attenuation by dust has as a function of wavelength. But even after proper 
extinction corrections are made, additional factors may lead to differences in the star forma- 
tion rate, including old stellar populations contributing to the infrared emission, incorrect 
Initial Mass Function (IMF) assumptions, differences in the relative stellar /nebular /dust 
geometries, contamination from Active Galactic Nuclei (AGN), and the differenti al impacts 
of metallicity and the ionization state o f the gas on the y a rious SFR indicators ( Kennicutt 



1998 



Xu et all |2003[ Ichoi et al.l booei ISalim et al.l boO?! : iLara-Lopez et al.l |2009|). While 



many of these discrepancies are due to the uncertainties in the conversion factors from lu- 
minosities to star formation rates, an additional possibility for discrepa ncy exists within the 
different survey selection criteria (ITresse et al.ll2002l : iBuat et al.ll2007bl ). More 'robust' star 
formation rate measures have rece ntly been devised using a combination of infrared -|- ultra- 



violet data or infrared -|- Ha data (iHirashita et al.ll2003t ICalzetti et al.ll2007l : iKennicutt et al. 



20071 : iPrescott et al.ll2007l : Kennicutt et al.ll2009t and references therein) because the infrared 



emission traces extinguished starlight while the ultraviolet and Ha emission trace relatively 
unextinguished starlight. 

These SFR indicators have bee n used in a large number of surveys to map the cosmic 



2006: Dale et al. 2008; Prescott et al. 



SFR at a variety of redshifts (e.g., ICowie et al.l 120041 : iLe Floc'h et al.l l2005l : iDoherty et al. 



20091) • The cumulative result of these efforts show that 



the cosmic SFRD has significantly evolved over the last sever al gigayears, perhap s by as much 
as a factor of 10 between ^ ~ 1 and the present epoch (e.g.. IWatson et al.ll2009l ). Therefore, 
it is not surprising that the galaxies we observe locally do not necessarily resemble those at 
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(Pei et al. 


1999; 


Davies et al. 


2009) 



ies (LIRGs) and Ultra-Luminous Infrared Galaxies (ULIRGs) are increasingly important as 
we look back to redshifts of z ~ 1 — 2 as bo th the IR and UV luminosity functions evolve 



toward a predorainance of LIRGs a t higher z (jPapovich et al.ll2006l : iReddy et al.ll2006l 12008 



Buat et al. 


2009; 


Salim et al. 


2009) 



the galaxy population as we move toward higher redshifts. Does this effect lead to a change 
in the average amount of obscuration by dust over time? Does the build up of metal content 



withi n galaxies as they age strongly affect their extinction by dust (e.g., iMannucci et al. 



20091 )? Not only do galaxies evolve with redshift, they each evolve on their own timescales 
(i.e., they run out of gas and cease to form stars, they are involved in mergers with other 
galaxies, etc.). These changes imply that the light that is detected from a given galaxy will 
be a function not only of its redshift, but also a function of parameters such as metallicity. 



star formation history, merger s 
of any supermassive black hole ( 


fcage (if applicable), and the relative h 


iminosity dominance 


Kennicutt 19981; ICalzetti & Heckman 


1999; 


Pei et al. 1999: 


Jansen et al. 2001: Dale et al. 2007; Wu et al. 2007; Eneelbracht et al 


2001 


^;and references 


therein) . 



Not only do SFR indicators tell us how many stars are being formed, but they can 
provide other physical insight when two different SFR indicators are coupled together. For 
example, the infrared-to-ultraviolet ratio is typically used to measure the amount of dust 
attenuation at ultraviolet wavelengths. Studies couphng infrared and ultraviolet data have 
shown that the slope of th e ultraviolet continuum is a useful prob e of the attenuation by 
dust in starburst galaxies (ICalzetti et al.lll994j : iMeurer et al.lll999l ). so simply tracking the 
infrared-to-ultraviolet ratio with lookback time should prove fruitful. Recent efforts, how- 
ever, lead to a conf using set of diverge n t results. Eviden c e from (res t frame) ultrayiolet- 
selected surveys by iBuat et al.l (l2007bl ). iBurgarella et al.l ( 120071 ) . and iReddy et al.l (120061 ) 
ai z ~ 0, 1, and 2, respectively, shows that the average dust att enuation within gal axies 
decreases with increasing redshift. Similar results were reported by lBuat et al.l (l2007al ) and 
Xu et al.l (120071 ) for LIRGs. This is the expected result if the dust within galaxies slowly 
builds up as th e effec ts of generations of stellar lifecycles accumulate over time. However, 
Le Floc'h et al.l (120051 ) find that the comoving infrared energy density for an infrared-selected 
sample evolves as (1 + z)^'^ at < z < 1, whereas the u lt raviolet only evolves a s (1 -|- z)"^-^ 
over the same s pan in redshift (ISchiminovich et al.ll2005l ). iTakeuchi et al.l (120051 ). using the 
same dataset as 
shift. While the 



Le Floc'h et al. 



Le Floc'h et al. 



(2005), rep ort a similar result o yer th e same range in red- 
(120051 ) and ISchiminovich et al.l (120051 ) results are based on 
samples of galaxies that are likely very different, together, they could imply that the av- 
erage ultraviolet attenuation by dust within galaxies increases with lookback time. This 
result would be counterintuitive if dust attenuation were only governed by metal abundance. 
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However, the result is understandable if the attenuation by dust increases with S FR (e.g., 
Wang k HeckmanlllQQQ : iHeckman et aPllQQsl : ICalzettil[200ll : iHopkins et aPboOlbh and the 



cosmic SFRD increases with 



ookb a ck time. Finally , a t hird possible r esult is reported by 



Bell et all fl2005r ). IChoi et all fl2006r ). IXu et al.l (120071 ). and lZheng et al.l (120071 ) who find dust 



attenuation values at higher redshifts (0.6 < z < 0.8) consistent with those found in the 
current epoch. However, it should be noted that all of these studies utilize some sort of IR 
selection. While not all of the studies using an IR selection come to this same conclusion, 
this possible selection effect should be looked at in more depth. 



Buat et al.l ( l2007bl ) have shown that ultraviolet and infrared surveys of the nearby uni- 
verse do not necessarily sample the same galaxy populations, and thus the above discrep- 
ancies may be in part due to sample selection. Th e ultraviolet continuum trac es the star 
formation rate on a time scale of ~10^ years (e.g., Ilglesias-Paramo et al.l 120041 ), while the 



infrar ed continuum traces star formation rates on a time scale of ~10'' years (e.g.. iKennicutt 



19981 ). This, as well as issues such as feedback from star formation rate evolution, Lyman a 
absorption, the evolution of dust grains and the evolution of the average dust temperature 
can all affect the infrared and ultraviolet luminositie s. Invoking data from a lar ge, uniform 
Ha-selected survey could help to clear this confusion ( Ilglesias-Paramo et al.ll2004j). Ha emis- 



sion traces star formation rates on a time scale of ;<10^ years (e.g., Ilglesias-Paramo et al. 
20041 ). making it an ideal comparison to the infrared continuum. Focusing just on this simi- 
larity in timescales, and ignoring evolutionary processes such as the dispersal of dust grains 
into the ISM, the ratio of Ha to IR should remain relatively constant over time within a 
galaxy. Recent work from the Spitzer Infrared Nearby Galaxy Survey (SINGS) find that 
the Ha;-|-24/xm data have an almost one-to-one correlation with that of extinction corrected 
Paschen a, considered an ideal star formation rate indicator as it directly traces H II re- 
gions and lies at a relativ e ly long wavelength to help minimize the eff ects of obscuration by 
dust JCalzetti et al.ll2007t IKennicutt e7al]l2007l : IPrescott et aDl2007t see Fig. 11 in Calzetti 
et al. 2007). This Ha-|-24/xm-based star form ation rate indica tor can be compared to the 
more commonly used indicator, L(Ha) alone (iKennicuttI Il998l ) to form a diagnostic sensi- 
tive to the amount of attenuation by dust within galaxies, a multi-wavelength diagnostic 
where all the inpu ts conveniently t r ace st ar formation on the same time scale. Although 
recent results from IKennicutt et al.l (|2009[ ) indicate that it can be just as effective to use a 
combination of the total infrared (TIR) emission -|- Ha or 8/im + Ha, TIR data are not 
frequently available for surveys of distant galaxies and 8/im data may be more problematic 



due to the variable nature of p olycyclic aromatic 



field hardness and metallicity (lEngelbracht et al. 



200a iRosenberg et al.ll2006l : IWu et al. 



2ooa 



lydrocarbon (PAH) emiss i on with radiation 



2005 



Dale et al. 



l2009al ) 



Madden et al. 2006; O'Halloran et al. 



The Wyoming Survey for Ha (WySH) is a large-area, ground-based, narrowband optical 
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imaging survey for Ha at redshifts z ~ 0.16, 0.24, 0.32, and 0.40 ( iDale et al.ll2008l ). In 



addition, this effort is being supplemented by the NEWFIRM Narrowband Ha Survey, a 



collab orative near-infrared imaging survey of Ha at redshifts of z ~ 0.81 and 2.2 (ILee et al. 



20091 : Lee et al., in prep; Ly et al., in prep.). Both surveys target key extragalactic deep fields 
where an abundance of multiwavelength imaging and spectroscopy are already available. In 
this work we utilize (Haobs + 24/im)/Haobs to directly gauge the amount of attenuation by 
dust in a large sample of galaxies that span redshifts out to 2; ~ 0.4. Our main goal is to 
robustly measure any changes in the average dust attenuation over cosmic time. Secondary 
goals include comparing the relatively new technique for assessing the attenuation by dust, 
by combining Ha and 24/im, with schemes based on the coupling of IR and UV. 

In this paper, we outline the datasets used in present the analysis and initial results 
from z ~ 0.16, 0.24, 0.32 and 0.40 in §31 and summarize our findings in §11 The cosmology 
used throughout this paper is Hq = 70 hjo km s^^ Mpc^^, Qm = 0.3, and Qx = 0.7. 



2. DATA 
2.1. Fields 

Results for the ELAIS-Nl (a ~ 16''11"^, S ~ +55°00' ) and Lockman Hole (a ~ lOMS'", 
S ~ +58°00' ) regions are presented here. These regions are popular targets for deep extra- 
galactic surveys since they are comparatively free of contamination by zodiacal and Galactic 
dust (/,y(100/im)~l MJy sr~^), as well as by bright stars and nearby galaxies. Two relevant 
surveys carried out in ELAIS-Nl and the Lockman Hole are the Spitzer Wide- Area Infrared 
Extragalactic (SWIRE) Survey (see § I2.2.2p and the Galaxy Evolution Explorer (GALEX) 
Deep Imaging Survey (see § I2.2.3p . Some thirty separate 18' fields were targeted by the 
WySH program in these two fields (see Figures [T] and [2]). Table [T] summarizes the total 
coverage at each epoch for both fields. 



2.2. Datasets 
2.2.1. Ha 

WySH is a ~4 square degree, ground-based, narrowband imaging survey for Ha emitting 
galaxies over the latter half of the age of the universe. WySH uses a pair of wavelength- 
adjacent narrowband filters at each redshift epoch (for z < 0.40) for improved contin- 
uum subtraction, relative to using a broadband filter to subtract the continuum. The 



Fig. 1.— The 24/im image from SWIRE of the ELAIS-Nl region {a ~ 16^11"^, 6 ~ +55°00' ) 
overlayed with circles to represent GALEX fields and squares to indicate the WySH fields. 
North is up and east is to the left. Dimensions of the SWIRE image are approximately 3° by 
3°, the GALEX overlays have diameters of approximately 1°, and the WySH overlays have 
dimensions of 18' by 18' . 
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Fig. 2. — The 24/im image from SWIRE of the Lockman Hole region {a ~ 10*^45™, 6 ~ 
+58°00' ) overlayed with circles to represent GALEX fields and squares to indicate the 
WySH fields. North is up and east is to the left. Dimensions of the SWIRE image are 
approximately 3.6° by 3°, the GALEX overlays have diameters of approximately 1°, and the 
WySH overlays have dimensions of 18' by 18' . 



Table 1. WySH Optical Imaging & Source Counts. 



z ELAIS-Nl (□") Lockman Hole (□") NB int. time (s) # Ha* # Ha^- + 24/imt> # Ha^- + 24/Ltmt' + UV^ 



0.16 1.38 1.01 1200 94 43 39 

0.24 1.35 0.83 3600 239 69 62 

0.32 1.33 0.87 6000 250 82 67 

0.40 052 OOO 9600 49 20 16_ 

Totals: ~ ~ ~ 632 214 184 



^Detections from WySH arc above 3o". 
''Detections from SWIRE are above 3<7. 
''Detections from GALEX are above 5(t. 
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epochs being observed at the Wyoming Infrar ed Observatory (WIRO) 2.3m telescope are 
z ~ 0.16,0.24,0.32, and 0.40 (jPale et aLll2008l ). This optical survey is being supplemented 
with collaborative near-infrared work at the Kitt Peak National Observatory 4m telescope 
at z ~ 0.81 and 2.2 by the NEWFIRM Narrowband Ha Survejil], which is observing through 
low airglow windows at 1.19 and 2.10/im (Lee et al., in prep.; Ly et al., in prep.). These data 
are being used to constrain the cosmic star formation rate density out to 2; ~ 2.2. WySH 
data are sensit ive down to a 3(7 survey depth of ^IM© yr~^ for z ~ 0.16,0.24,0.32 and 
0.40, using the iKennicuttI (119981 ) equation to convert from L(Hq;) to SFR(Hq;). This value 
has not been corrected for dust attenuation. For a full description of the optical survey and 
pre liminary re s ults for the SFRD and luminosity function at z ~ 0.16, 0.24, 0.32 and 0.40, 



see 



Dale et all (120081 . l2010[ ). 



2.2.2. Infrared 



SWIRE (ILonsdale et al.ll2003l ). one of the Spitzer Legacy Projects, was designed to 
image ~49 square degrees of the sky, as well as detect very distant galaxies (out to z ~ 3) 
in order to look at different properties of galaxy (and star) formation and evolutior^. SWIRE 
uses all seven imaging bands o n Spitzer (3.6, 4.5, 5. 8 and 8.0/im with IRAC and 24, 70 and 
160yum with MIPS). Using the ICalzetti et all (|2003) relation between L^(24/im) and SFR, 



SFR(Mq yr~i) = 1.27 x 10"3^[z/L^(24/im)(erg s'^)] 



0.8850 



we find that, for a 3a 24/im survey depth of 135/iJy, at z ^ 0.16,0.24,0.32 and 0.40, 
galaxies with SFRs of approximately 0.22, 0.50, 0.86 and 1.4 Mq yr~^ (respectively) can be 
detected. The SWIRE survey covers 9.2 square degrees in the ELAIS-Nl field and 11 square 
degrees in the Lockman Hole field (see Figures [1] and [2]) . 



2.2.3. Ultraviolet 



GALEX (IMartin et al.ll2005l ) is a space-based mission designed to study the ultraviolet 
properties of galaxies in the local univers^. Utilizing both the far- and near-ultraviolet 



^ http: / /newfirm. ociw.edu/wiki/Main_P age 
2www.ipac.caltech.edu/SWIRE/ 
www. galex.caltech.edu 
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(FUV, NUV) detectors (1538.6A and 2315. 7A, and FWHMs of 269A and 616A respec- 
tively), the Deep Imaging Survey pointings have exposure times of 30ks, as opposed to the 
O.Iks exposure times per pointing for the GALEX All Sky Survey. Deep Imaging Survey 
pointings include ELAIS-Nl and the Lockman Hole. The Deep Imaging Survey has 5a 
limiting magnitudes (in m^^) of mAr[7y=24.4 and mFuv='24:.8, angular resolutions of 5.3" 
("NUV) and 4 . 2" (FUV) and samples gala xies with an average redshift of < z >= 0.85 
( Martin et al.l 2005 : Morrissey et al. 2007 ). Using the FUV limiting magnitude and the 
KennicuttI ( IQQsi l relation between UV luminosity and SFR, we find that, for a 5cr survey 
depth of = 4.365 x lO'^^ W m'^ Hz-\ at z ~ 0.16,0.24,0.32 and 0.40, SFRs of 
approximately 0.04, 0.11, 0.20 and 0.35 Mq yr~^ (respectively) can be detected (assum- 
ing no attenuation by dust). It is also important to note that the GALEX Deep Imag- 
ing Survey is confusion lim ited and that standard pipeline results are not reliable for fain t 
i^Nuv > 23 mag) sources (jBurgarella et al.ll2007l : iMartin et al.ll2007l : IZamojski et al.ll2007l ). 
This may slightly underestimate dust attenuation using the i'L^{24:fj,m) / uL^^FUV) indica- 
tor for the UV faint galaxies. These UV faint galaxies follow the same trend as the brighter 
{rriNuv < 23 mag) UV galaxies, and make up ~8% of the entire sample ([0, 5, 8, 2] sources 
at z ~[0.16, 0.24, 0.32, 0.40]). For the joint WySH/GALEX Deep Imaging Survey analysis 
carried out here, a total of 12 ~1 degree fields covering the ELAIS-Nl region and 18 ~1 
degree fields covering the Lockman Hole region (see Figures [1] and [2]) are used. 



2.2.4. Merged Catalog 

The Ha data from WySH were used for the initial sample selection, providing a sample 
of 94 + 239 + 250 + 49 galaxies at 2 ~ 0.16,0.24,0.32, and 0.40, respectively based on 
a 3cr detection limit. This data was first merged with a separate 24/im SWIRE catalog 
with a depth of 3a. Then, using the GATOR toofl (provided through the NASA/IPAC 
Infrared Science Archive (IRSA)), the Spring 2005 SWIRE Spitzer catalogs for both the 
ELAIS-Nl and the Lockman Hole fields were accessed. These catalogs contain data in the 
four IRAC bands, as well as the MIPS 24/im band. All sources in the catalog contain 3.6/im 
and 4.5/im data, but not all sources were detected at the necessary level in the remaining 
bands. In GATOR, we uploaded the WySH source coordinates separately for each epoch 
{z ~ 0.16, 0.24, 0.32, 0.40), and searched for coordinate matches with SWIRE sources within 
a 3" radius to create a WySH + SWIRE database. To create the WySH + SWIRE + GALEX 
merged catalog of sources, the WySH -|- SWIRE catalogs were coordinate matched with the 
GALEX catalogs from GALEX Release 5 (GR5) using a search radius of 3". The distance 



http://irsa.ipac. caltech.edu/applications/Gator/ 
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between any two Ha sources is larger than the point-spread functions (PSFs) of both Spitzer 
and GALEX, so the ~6" resolution of MIPS24 and GALEX imaging usually does not result 
in multiple galaxies contributing to a single IR or UV flux. However, if multiple SWIRE 
or GALEX counterparts are found, the closest one is used. For Ha sources having Ha 
neighbors within 3" (7 sources), there are 3 UV sources which appear as extended (blended) 
sources, which accounts for < 2% of the entire sample. The sample is further reduced in 
size by requiring sources to have data in both UV bands as well as the 24/im band. Both 
the WySH + SWIRE and WySH + SWIRE + GALEX catalogs are used in this research 
where appropriate. Ultimately, this can be considered a 24/xm selection. The sample sizes 
are listed in Table [H 



2.3. Corrections 



2.3.1. k Corrections 



For k corrections, the central redshift value for each epoch is used rather than the redshift 
for each source within a given epoch since Sz per epoch is ~0.02 a nd the redshifts for this 



sample have been determined photometrically (IRowan- Robinson et al.ll2008l ). The ultraviolet 
data for the sample were k corrected by calculating the slope between the observed FUV and 
NUV fluxes and by extrapolating from that slope to the rest frame FUV and NUV fluxes. 
The FUV k corrections are modest and have median values of fy{FUVy'"^'^/fy{FUV)°^'^ = 
[0.91,0.89,0.86,0.85] at ^ ~ [0.16, 0.24, 0.3 2, 0.40]. Comparing this to local data fro m the 
Spitzer Local Volume Legacy (LVL) survey (IKennicutt et al.ll2008l : iDale et al.ll2009bl :Lee et 
al. 20 09, in prep.; for Ha, 24/ xm and UV data, respectively), and the Nearby Galaxy Survey 
(NGS: lGil de Paz et al.ll2007l ) in Figure[3l the ultraviolet continuum of late-type galaxies (T- 
type > 3) is also relatively flat (median values along the y-axis are 0.06 dex for the Spitze r 
Local Volume Legacy s urvey and 0.02 d e x for NGS). (Also see Figure 2 of iDale et al.l (120071 ) 
and Figures 5 and 6 of iGil de Paz et al.l (120071 ).) For the ii ifrared data, the k co rrections are 
derived based on the average star- forming SED models of iDale fc Heloul (|2002| ). The 24/im 
A; corrections are /^(24)*™7/^(24)°''^ = [1.57,1.90,2.25,2.63] at ^ ~ [0.16, 0.24, 0.32, 0.40]. 



2.3.2. AGN Corrections 



The empirical method of iLacy et al.l (120041 ) is used to identify possible AGN w ithin 
our sample (see Figure Hj). The different "sequences" identified in iLacy et al.l (120041 ) are 
not immediately apparent in Figure H] since the large Sloan Digital Sky Survey (SDSS) 
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Fig. 3.— The ratio of \og{uF^{FUV) /uF^{NUV)) v. \og{uL^{NUV)). Small pluses repre- 
sent data from late type (T-type > 3) galaxies in the Spitzer Local Volume Legacy survey and 
NGS, while small diamonds represent early type (T-type < 3) galaxies in the Spitzer Local 
Volume Legacy survey and NGS. These local data show that log^uFj^^FUV) /uF^^NUV)) 
remains relatively flat with median values (for the late type galaxies) of 0.06 dex for the 
5'p«t2;er Local Volume Legacy survey and 0.02 dex for NGS, with respect to logi^uL^^NUV)) . 
Rest-frame data from WySH are also shown. Triangles represent data at z ^ 0.16, squares 
represent data at z ~ 0.24, crosses represent data at z ~ 0.32 and diamonds represent data 
at z ~ 0.40. 
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dataset is not also included here as in Figure 1 of iLacy et al.l (|2004[ ). iLacy et al.l (120041 ) 
describe one "sequence" to have blue colors in Ly(5.8/im)/Ly(3.6/im) and very red col- 
ors in Li,(8.0/im)/Ly(4.5yum), which they describe as likely to be low-redshift {z < 0.2) 
galaxies whose PAH emission bands have not been significantly redshifted out of the 8.0/im 
filter. Within this same sequence but at higher redshifts, the 7.7/im PAH feature is no 
longer observed with the S.O/im filter, which would give those sources a bluer color in 
Li,(8.0/im)/Lj,(4.5/im), and a redder color in Lj,(5.8/im)/Ly(3.6/im) since the 3.3/im PAH 
feature will have redshifted out of the 3.6nm band . In F igure [H it is clear that most of the 
sources in this sequence as defined by iLacy et al.l (|2004( ) are found at either z ~ 0.16 or 
z ~ 0.24, however, the sources at z ~ 0.32 and z ~ 0.40 start to span the gap between the 
two sequences due to the e ffects of redshift a nd lie closer to, or within, the enclosed region. 
The other "sequence" that iLacy et al.l (|2004[ ) describe has red colors in both filter pairs; this 
is where they find SDSS- and radio-s elected qua s ars. T he dashed lines in Figure S] enclose 
this region , as de fined empirically by iLacy et al.l ( 120041 ). Regarding the WySH sample, the 



Lacy et al.l ( l2004l ) empirical classification of AGN is used to conservatively identify the three 



right-most sources within the box in Figure H] as AGN. These three sources are removed 
from the remainder of this study. The identities are less clear for the other sources near the 
AGN/star- forming boundaries, which are indicated by dashed lines. It is possible that some 
of these sources are composite sources (i.e., a combination of AGN and star formation within 
the galaxy). Other sources, particularly at 2; ~ 0.32 and z ~ 0.40, may have simply been 
redshifted into this region based on the previous description of the "sequences". Without 
spectra, it is impossible to confirm whether or not the sources in question are AGN, star 
forming or a composite, and these sources near the AGN/SF boundaries are not removed 
from the sample. 



3. RESULTS 

3.1. Dust Attenuation via 24^m & FUV 

Very few sources in this study have detections of far-infrared emission, so the traditional 
L(TIR)/L{FUV) tracer of the attenuation by dust cannot be used, and is replaced here with 
i'Li,{24fim)/i'L^{FUV) as proxy. The rest-frame 24/im emission is approximately 10% of 
the TIR emission on average for normal star forming galaxies, but the dispersion is large. 
For normal star-forr ning galaxies at the epochs studied here, the full range in 24/im/TIR 



spans a factor of 2 ( jPale et al.ll2005t see Figure 16). Figure |5] shows the dust attenuation 



diagnostic, i'L^{24fim)/i>Li,{FUV) as a function of the Ha+24yum star formation rate from 
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Fig. 4. — An IRAC color-color plot of the WySH sources. Triangles represent sources at 
z ~ 0.16, squares represent sources at 2; ~ 0.24, crosses represent sources at z ~ 0.32 
and diamonds represent sources at 2; ~ 0.40. The region enclosed by the dashed hne shows 
the empirically de termined color criteria used to distinguish possible AGN within the sample 
jLacv et al.ll2004[ ). 
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Kennicutt et al.l f l2009h : 

SFR{MqYt~^) = 7.9 X 10-^^[aL{Raobs) + bL^{2Afim)]. 



(2) 



where both lu minosities are in erg s~ \ and L^{2Afim) is expressed as uL^,. The best fit to 
the data from iKennicutt et aLl (120091 ) gives b/a = 0.020, which is based on the integrated 
fluxes of a large sample of galaxies. Figure [5] follows the trend that galaxies with larger star 



format ion rates tend to be more dusty, as inferred from the 24/im/FUV ratio ( iHopkins et al. 



2001al ). Looking at only the star formation rate portrayed along the x-axis, a trend toward 
higher star formati on rates at higher r edshifts is seen, which has been shown by a number of 



2008: Prescott et al. 



other studies fe.g..lCowie et al.l 120041 : iLe Floc'h et al.l 120051 : jPohertv et al.ll2006l : jPale et al. 



2009! 



|. This evolution could partially be due to sample selection whereby 



only the brighter galaxy subsets are observed at higher redshifts. The data from all four 
epochs appear to have roughly the same average attenuation by dust, as indicated by the 
large symbols in Figu re |5] (and th e infor mation portrayed in Figure[7]). This result agrees with 
the recent results of iBuat et al.l ( l2009l ). While the ensemble average is relatively constant 
from epoch to epoch, each epoch separately exhibits a larger attenuation by dust for higher 
star formation rates. Hence, given that the average galaxy star formation rate increases with 
lookback time, an epoch to epoch comparison at a flxed star formation rate suggests a mild 
decrease in dust attenuation with redshift. The argument of a decrea se in the amount o f 
attenuation by dust with lo okback time is cons is tent with re s ults fro m lBuat et al.l (l2007a ) 



and c ombined results from iReddy et al.l (120061 ): iBuat et al.l (l2007bl ) and iBurgarella et al. 
(120071 ) ■ who have found dust attenuation values decrease from z ~ to z ~ 2, as discussed 
previously. 



3.2. Dust Attenuation via Ha & 24yum 



The empirical argument for a decrease in dust attenuation is perhaps more readily 
apparent from the data portrayed in Figure El Figure [HI shows the dust attenuation di- 
agnostic, [aLfH a ohs) + &Lj,(24yum)]/L(Ha;obs) versus the star formation rate indicator of 



Kennicutt et al.l (120091 ) (Equation |2]). Here, again, the expecte d trend that galaxies with 
higher star formation rates tend to be more dusty can be seen (IHopkins et al.ll2001al ). al- 
though, in this case the trend is seen most strongly for a given redshift rather than for the 
entire dataset. The right hand ax:is tra nslates this diagnostic into a value of magnitudes of 
Alia extinction (IKennicutt et al.l l2007l ) . The large symbols in Figure E] represent the aver- 
age values for each epoch. As in Figure |5], the average dust attenuation remains relatively 
constant from epoch to epoch. Moreover, it is likewise more illuminating to compare the 
attenuation from epoch to epoch at a flxed star formation rate. For a given star formation 
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Fig. 5. — A proxy for the traditional measure of dust attenuation, uL ,, (24Mni) / i>L,y (FUV) (i n 
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Triangles represent sources at 2; ~ 0.16, squares represent sources at 2 ~ 0.24, crosses 
represent sources at z ~ 0.32 and diamonds represent sources at 2; ~ 0.40. Large 
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rate, the typical amount of attenuation by du st decreases with inc reasin g redshift, echo- 



i ng th e combined results of iReddy et al.l (120061 ) . iBuat et al.l (l2007bl ). and iBurgarella et al. 
( 120071 ). While these studies use a UV selection, no compara ble studies exist us ii ig an Ha 



selection. F or work based on ot her selections, see, for example, iBuat et al.l (l2007al ). IXu et al. 



( 120071 ) , and IZheng et al.l (120071 ) . These comparisons do not change the result that while the 
ensemble average of dust attenuation does not change as a function of the redshifts sam- 
pled here, we expect higher redshift galaxies to show higher values of attenuation since they 
are typically exhibiting larger star formation rates. However, we find that higher redshift 
galaxies show lower values of attenuation at a given luminosity (i.e., star formation rate). 



3.3. Comparison of Dust Attenuation Diagnostics 



To compare Figures [5] and [6], the correlation coefficients for each diagnostic were com- 
puted. The \og{h'L,y{24:fim) /uL^i^FUV)) diagnostic gives r = 0.545 and the log([aL(Haobs) + 
6L,^(24yum)]/L(Haobs)) diagnostic gives r = 0.452. To investigate the stati stical significance 
of the two diagnostics, the program R§| ( R Development Core Team 2009 ) was used to an- 



alyze the data with a linear model. Based on the multiple statistic, ~0.11% of the 
variability in log(z/Lj^(24/im)/z/L,y(Ff/V^)) can be explained by the star formation rate, and 
~0.01% of the variability in log([aL(HQ;obs) + &-^i/(24/im)]/L(Haobs)) can be explained by 
a linear association with \og{SFR). It can be argued that the latter diagnostic is slightly 
more correlated with SFR because the [aL(Haobs) + bLi^{24:fim)] parameter appears on both 
axes, however, based on the statistical analysis described above, this positive correlation is 
not seen in the data. 

The dust attenuation seems to be higher when looking at the line emission from Ha 
(Fig. [6]) versus the UV continuum (Fig. [5]). However, it is more informative to look at the 
amount of attenuation by dust determined from both diagnostics on the same scale {Av^„„t 
(continuum) in magnitudes) with respect to redshift (see Figure [7]). In order to convert 
pLi, i24nm) /uLi, ( FUV) to Av^^^^-, z/Lj,(24/im) was multiplied by 10 to approximate TIR 
(jPale et al.ll2005l :Fig. 16). (This is a large so urce of uncertain ty, approximately a 
factor of two.) The following prescription from lBuat et al.l ( l2007al ) was used to determine 
^Vcont (ill magnitudes) from TIR/ FUV: 



A{FUV) = -0.0333 log 



'FUV 



+ 0.3522 log 



L 



TIR 
'FUV 



+ 1.1960 log 



■'TIR 
'FUV 



+ 0.4967, 
(3) 



http://www.r-project.org/ 
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where the numerator is t aken from Kennicutt et al. ( 20091 ) as a function o f the HQ;+24/zm 
star formation rate from Kennicutt et al. ( 20091 ). and h/a = 0.020. ( Kennicutt et al. 



20071 ) is plotted on the right hand axis. Symbols are the same as in Figure [51 
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A{FUV) 



(4) 



The reddening curve from iLi fc Draind (120011 ) converts from A^^ to Ay^^ (emission line) for 



the [aL(Haobs) + &i^,y(24/im)]/L('Hanhs ) diagnostic. The value R'y = 4.05 is used to convert 



from Av to A 



1.8Av,,^t (ICalzetti et al.ll2000l : ICalzettill200l[ ). As can be seen in 



Fig. [71 when the two diagnostics are placed on the same Av^^„t scale, they agree within the 
errors. This consistency holds whether the errors are determined through population statis- 
tics or formal measurement uncertainties. Both diagnostics show no noticeable evolution in 
the amount of dust attenuation as a function of redshift. 

There are additional challenges to comparing uLy{24:iJi,m) / vLy{FUV) and [aL(Ha;obs) + 
6L,y(24/im)]/L(Haobs)- The Ha and the infrared continuum emission both trace star forma- 
tion on time s cales of < l{fYi, while the ultraviolet cont i nuum traces SFRs on time s cales of 
~ 10 V (e.g., lKennicuttiri998[ iMevnet k Maederl I2OO0I : lldesias-Paramo et all 12004 ). Thus 
the [aL(Haobs) + &-^^i/(24/im)]/L(Haobs) is naturally a 'cleaner' diagnostic of attenuation by 
dust since all factors involved arise from similar regions, which are physically separated 
from much of the UV emission , as well as similar time scales (for a counter- argument, see 
Gordon et al.l I2OOOI : iBuat et al.l 120051 : ICortese k. Hughes! |2009[ ) . An additional challenge to 
interpreting global i'Ly{2'i^m) / vL^i^FUV) measures lies in the average dust temperature of 
a galaxy. The 24/im emission is not only sensitive to the amount of dust, which directly 
leads to attenuation, but it is also sensitive to the temperature of the dust. In the dust 
models of iDraine fc Lil ( 120071 ). the flux at 24/im largely represents emission from very small 
grains with effective radii of 15-40A. These intermediate-sized dust grains are responsible for 
a widely varying portion of the bolometric infrared luminosity depending on the intensity 
of the interstellar radiation field. As the intensity of the radiation field increases linearly, 
the 24/im emission from the dust grains will increase more r apidly, appearing ov er-luminous 
with respect to the linear increase of the radiation field (see ICalzetti et al.l 120071 ). Naturally, 
any monochromatic indicator of the dust emission is limited in its ability to effectively en- 
capsulate the emission processes of all dust grains that emit between 1 and lOOO/xm, and the 
24/im bandpass is no exception. 



3.4. Comparison with Dust Attenuation Results at z 







It is also important to compare the WySH dataset to a dataset at 2; ~ 0, in order 
to increase the redshift baseline to see how the amount of attenuation by dust evolves from 
z ~ to 2 ~ 0.40. Using the same diagnostics as in Figures [5] and El data from the Spitzer 



Loca l Volume Legacy s urvey and SINGS are added to Figures M and M ( IKennicutt et al. 



2008uDale et al.ll2009bl :Lee et al. 2009, in prep.; for Ha, 24/im and UV data, respectively). 
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These samples were chosen because of their multi-wavelength datasets that somewhat mirror 
the data already being used in this study. In particular, the Spitzer Local Volume Legacy 
survey also contains a l arger number of lo w-luminosity galaxies than would be detected at 
the redshifts of WySH (jPale et al.ll2009bl ). It is clear from Figures |H] and [H] that the local 
population of galaxies provides a natural extension to lower SFRs and lower levels of dust 
attenuation. Unfortunately, where the WySH and local populations overlap in SFR, there 
are not enough data to infer any meaningful conclusions. The emphasis of this comparison 
is on the fact that the dust attenuation does, in fact, increase with SFR. 



4. SUMMARY 



The study of the evolution of dust attenuation within galaxies over cosmic time has led 
to a set of divergent resu lts. Indeed, the at tenua tion by dust has be e n found to decreas e 
from z ~ to ^ ^ 2 in lSuat et all J2007ah and iReddv et all J2006h: lBuat et all J2007bh: 
Burgarella et al.l (120071) (combined), to increase from z ^ to z ^ 1 



2005 



Le Floc'h et al. 



from z ~ 0.6 — .8 to the present epoch (IBell et al. 



Schiminovich et al.l l2005l : iTakeuchi et al.l 12005 ) and to remain essent i ally unchanged 



2005; Choi et al. 2006; Xu et al. 2007 



Zheng et al.ll2007l ). These disparate results could be the effect of different sample selection 
critiera, with one study focusing on UV-bright galaxies, another on ULIRGS, etc. The 
selection criteria for the various studies are summarized in Table |2J There does not appear to 
be any one single selection criteria that maps to any of the three above-mentioned categories 
of results (dust attenuation increases with z, decreases with z, or remains unchanged), so it 
is difficult to determine whether or not (or to what extent) the selection criteria drive the 
conclusions made in the previous efforts. 
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Table 2. Summary of previous studies of dust attenuation as a function of increasing 

redshift. 



Method^ Dust Attenuation Amplitude'' Source Selection Reference 



dust attenuation decreases 



LD 


< ^PUV >~ 3 mag and decreases 


24/im selection (z < 1); LIRGs only 


Buat et al. f 2007a) 




from 2: = to 2 = 0.7 






LF 


j4fuv decreases from ~ 3.4 to 1.7 


mag separate FUV & FIR selections (z ~ 0) 


Buat et al. ('2007b) 




from z = to z = 2"^ 






LF 


Apuv decreases from ~ 3.4 to 1.7 


mag UV selection of Lyman break galaxies 


Burgarella et al. (2007) 




from z = to z = 2"^ 


(z ~ 1) 




LF 


j4puv decreases from ~ 3.4 to 1.7 


mag BM/BX (optical) criteria similar to Lyman 


Reddv et al. f2006) 




from z = to z = 2"^ 


break galaxies atz ~ 3(z ~ 2) 




Ind 


decreases 1.5 mag 


24/im selection with Ha & UV detections 


this study 




from z ~ 0.16 to z ~ 0.40 


required 




dust attenuation increases 


LF 


co-moving IR energy density 


24/im selection (z < 1) 


Le Floc'h et al. ("2005) 




{l + z)3-9±0'4 (0<z<l)'i 






LD 


Pl500 = (l + ^)2'5±°-^toz~ld 


UV selection (z < 3) 


Schiminovich et al. f200E 


LD 


Pdust / Pfuv increases from ~ 4 


FUV & IR luminosity functions from 


Takeuchi et al. (2005) 




to ~ 15 (z = to z ~ 1) 


previous works 




dust attenuation same as at 2 ~ 


Ind 




NIR + MIR selection (2 ~ 0.8) 


Choi et al. C2006) 


LD 




24/im selection (z ~ 0.7) 


Bell et al. (2005) 


Ind 




separate NUV & 24/im selections of LIRGS 


Xu et al. (2007) 






& ULIRGS (z - 0.6) 




LD 




COMBO-17 photo-z (0.6 < z < 0.8) 


Zhene et al. (2007) 



''Results are derived from either a luminosity function (LF), a luminosity density (LD) or a collection of individual galaxies 
(Ind). 

''Amplitude of dust attenuation evolution & what range the evolution is measured over. 



'^The resu lts fromlBuat et alj ll2007lJ l : liurgarella et"al] l|2007b and lReddv et alj ll2006h are combined to determine the results 
presented in lBurgarella et al. out to z ~ 2. 

'^The results from iLe Floc'h et al. and ISchiminovich et al. I ll2005l'l are combined to reach the conclusion from 

iLe Floc'h et al.l l l2005h that dust attenuation increases at higher z. 



to 
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In an effort to decipher this problem, we present multi-wavelength results that explore 
the evolution of the amount of dust attenuation from z ~ 0.16 to z ~ 0.40. We use 
Ha data from WySH for a uniform initial source selection, and further constrain our sample 
by selecting according to detections available in IR data from SWIRE and UV data from 
GALEX. Using a traditional measure of dust attenuation, i/L,^(24/xm)/i/L,^(Ft/y), we show 
that over this range in redshift, there is no noticeable evolution in the ensemble average 
amount of attenuation by dust of our sources. Looking over this same range in redshift, 
with a more recent measure of dust attenuation, [aL(Haobs) + &-^i/(24/im)]/L(Haobs), "we also 
show that there is no noticeable evolution in the ensemble average amount of attenuation by 
dust. However, looking epoch by epoch for a given star formation rate, a mild decrease in 
dust attenuation with lookback time can be seen. Adding data from a couple of local galaxy 
samples, spanning both a volume-complete sample {Spitzer Local Volume Legacy survey) 
and one that is more geared toward IR-bright galaxies (SINGS) shows no (or very little) 
evolution in dust attenuation in galaxies from z ~ to ^ ~ 0.40. Two general trends 
found in the data are that the dust attenuation increases with star formation rate and that 
galaxies were probably intrinsically more luminous at times in the past. 

Future efforts will incorporate Ha data at 2; ~ 0.8 and 2.2 from the NEWFIRM 
Narrowband Ha Survey. We also plan to extract our own flux measurements from the 
GALEX data. The NEWFIRM Narrowband Ha Survey collaboration is also obtaining 
spectral data for the Ha sources that are identified at 2; ~ 0.8 and 2.2 in order to measure 
the dust attenuation independently using the Ha/H/3 ratio (Momcheva et al., in preparation). 
This will allow for a confirmation of the use of the empirical star formation rate equation 
determined locally with SINGS data (Eqn. |2]) as a dust attenuation diagnostic tool. 
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